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In this Letter we report the first results on % , K*, p and p production at mid-rapidity (\y\ < 0.5) in 
central Pb-Pb collisions at y / s r NN = 2.76 TeV, measured by the ALICE experiment at the LHC. The pj 
distributions and yields are compared to previous results at ^/snn = 200 GeV and expectations from 
hydrodynamic and thermal models. The spectral shapes indicate a strong increase of the radial flow 
velocity with y'sNN, which in hydrodynamic models is expected as a consequence of the increasing 
particle density. While the K/n ratio is in line with predictions from the thermal model, the p/n ratio 
is found to be lower by a factor of about 1.5. This deviation from thermal model expectations is still 
to be understood. 
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High-energy heavy-ion collisions offer the unique possibility to study nuclear matter under extreme con- 
ditions, in particular the deconnned phase (Quark-Gluon Plasma, QGP CD El IS), which has been pre- 
dicted by lattice QCD @. The transverse momentum (pj) distributions and yields of identified particles 
are instrumental to the study of the collective and thermal properties of this matter. Results from lower 
energies, in particular from the Relativistic Heavy-Ion Collider (RHIC, ^/^nn = 200 GeV), have shown 
that the bulk matter created in high-energy nuclear reactions can be quantitatively described in terms 
of hydrodynamic models. The initial hot and dense partonic matter rapidly expands and cools down, 
ultimately undergoing a transition to a hadron gas phase [5]. The observed particle abundances were 
described in terms of thermal models. Relative particle abundances in thermal and chemical equilibrium 
are governed mainly by two parameters, the chemical freeze-out temperature T c i, and the baryochemi- 
cal potential Hb, where the latter describes the net baryon content of the system ||6l|7l[H0. Measured 
particle yields in heavy-ion collisions at RHIC, as well as SPS and AGS, are consistent with equilib- 
rium populations, allowing the extraction of both model parameters from fits to the measured particle 
ratios |0 [71 [TUl. It has been argued that interactions modifying the relative abundances of particle 
species are negligible in the hadronic phase iTTTl PT21 , and that T c f, can be linked to the phase transition 
temperature |[L3"ll . Particle momentum distributions reflect the conditions later in the evolution, at the 
"kinetic freeze-out" from the hadron gas phase, when elastic interactions end Ifl4ll . The pj distributions 
encode information about the collective transverse expansion (radial flow) and the temperature T^m at 
the kinetic freeze-out [ T51 [T6ll . The collective expansion is driven by internal pressure gradients and 
quantified by the average transverse velocity (fh)- Based on the success of the thermal and hydrody- 
namic models and based on the trend of the model parameters with ^Ann, predictions were formulated 
for higher energy Q[T7l. With the advent of the LHC, it became important to check these expectations 
in the new energy regime. In this Letter, we present the first results on n, K, and p production in 0-5% 
central Pb-Pb collisions at -^snn = 2.76 TeV, measured by the ALICE experiment at the LHC. Previ- 
ous results in pp collisions are reported in |[T"8l . The measurement spans the pj range from ~0.1 up to 
~4.5 GeV/c at central rapidity (|j| < 0.5). 

The central tracking and particle identification (PID) detectors cover the pseudorapidity window 1 77 1 < 
0.9 and include, from the innermost outwards, the Inner Tracking System (ITS), the Time Projection 
Chamber (TPC), the Transition Radiation Detector (TRD) and the Time- Of -Flight array (TOF) ITT91I20T1 . 
The central detectors are operated in a 0.5 T solenoidal field. The moderate field, together with a low 
material budget (X/Xo « 0.1 for a track going through the TPC) permits the reconstruction of low pi 
tracks. The data sample consists of 4 M minimum bias events. A pair of forward scintillator hodoscopes, 
the VZERO detectors (2.8 < r\ < 5.1 and —3.7 < 17 < —1.7), was used for triggering and for centrality 
determination ET1 1221 |23l . The data were collected using a minimum bias trigger requiring a combi- 
nation of hits in the ITS and in the VZERO, a condition fully efficient for the event sample discussed 
here ||23l . The 0-5% most central collisions were selected using the signal amplitudes measured in the 
VZERO, fitted with a Glauber model [23 ]. Background events caused by beam-gas interactions, parasitic 
collisions and electromagnetic processes were rejected using timing cuts on the VZERO and two neutron 
Zero-Degree Calorimeters located on either side of the interaction point, at 1 14 m distance ll2Tll22ll23ll . 
The measurement was performed in three independent analyses, each one focusing on a sub-range of 
the total pj distribution, exploiting the capabilities of the individual detectors and specific techniques 
to optimize the signal extraction. The ITS is a six-layered silicon detector. The two innermost layers, 
based on silicon pixels, are also used in the online trigger as mentioned above. The four outer layers, 
consisting of drift and strip detectors, provide identification via the specific energy loss. Using the ITS 
as a standalone tracker enables reconstruction and identification of low-momentum particles that do not 
reach the TPC, in the pj ranges 0.1-0.5, 0.2-0.5, 0.3-0.5 GeV/c for %, K, and p. For each track, at least 
three dE/dx samples were required. Only the lowest two were used in a truncated mean procedure, 
leading to an ~10% resolution. The particle identity was assigned according to the distance from the 
expected energy loss curves, weighted by the resolution. This procedure results in asymmetric ranges 



2 



The ALICE Collaboration 



around the curves for %, K, and p, reflecting the Landau tails in the detector response, which are not fully 
suppressed by the truncated mean. An additional 2 a cut was applied in the case of pions, to remove the 
contamination from electrons. The residual misidentification (< 10% for kaons and negligible for pions 
and protons) is corrected using Monte Carlo. The other analyses combined the tracking information from 
the ITS, TPC and TRD ("global tracks") ||T9ll . The TPC identifies particles via the specific energy loss 
in the Ne-C02-N2 (85:10:5) gas mixture. Up to 159 samples are measured, but only the lowest 60% 
are used in the analysis. This truncated mean procedure leads to a Gaussian distribution with an ~6.5% 
resolution. The TOF is placed at a radius of 370-399 cm. It measures the time-of-flight of the particles, 
allowing hadron identification at higher pj. With a total time resolution of about 85 ps, PID is possible 
up to pj = 3 GeV/c for % and K, and 4.5 GeV/c for p. In the intermediate pj range (0.3-1.5 GeV/c, 
0.3-1.3 GeV/c, 0.5-2.4 GeV/c, for 71 and K, and p), the identification was based on the combined TPC 
and TOF signals. The ranges were chosen such that the contamination from misidentification of other 
species is negligible. It was required that the particles are within 3 a from the expected dE/dx and 
time-of-flight values. The TOF information was included starting at pj = 0.65, 0.6, 0.8 GeV/c for %, K, 
and p. In the third analysis, for pj > 0.5 GeV/c, the TOF signal alone was used for identification. The 
time-of-flight distribution was fitted for each pj bin with data-derived templates for %, K, and p, allowing 
to extract the particle yields when the separation is as low as ~2 a 11241 . This fit was repeated for each 
mass hypothesis, after applying the selection \y\ < 0.5 based on the mass assumption under study [24]. 

In this Letter, results for "primary" particles are presented, defined as prompt particles produced in 
the collision, including decay products, except those from weak decays of strange particles. Both ITS 
standalone and global tracks provide very good transverse impact parameter resolution relative to the 
primary vertex (DCA^), of order 200 at pj = 300 MeV/c and 35 at pj = 5 GeV/c, allowing 
to separate primary and secondary particles. The residual contamination was estimated from data by 
fitting the DCA xy distribution with three Monte Carlo templates: primary particles, secondaries from 
material and secondaries from weak decays 11251 [T8l . This contamination can reach 35% at pj = 300 
MeV/c for protons. It quickly decreases with increasing pj, becoming negligible at pi ~2.5 GeV/c. 
The efficiency correction and the templates used in the correction procedure were computed with 1 M 
Monte Carlo events, using the Hijing ll26l event generator (tuned to reproduce the dA^h/drj measured 
for central collisions ll22l ). transported through a GEANT3 IT271 model of the detector. The results of 
the three analyses were consistent in the regions of overlap and therefore combined using the (largely 
independent) systematic uncertainties as weights. 

The main sources of systematic uncertainties are summarized in Table [T] The uncertainties due to the 
secondary subtraction procedure were estimated for all analyses varying the range of the DCA rv fit, using 
different track selections (for instance using TPC-only tracks), applying different cuts on the longitudinal 
DCA Z and varying the composition of the Monte Carlo templates used in the fit. The uncertainty on the 
energy loss correction was estimated in a simulation with the material budget scaled by ±7%. In order 
to account for the uncertainties due to poorly-known hadronic interactions with the detector material, 
different transport codes (GEANT3, GEANT4 and FLUKA (29l) were tested. The interaction cross 
section used in the models for %, K, and p were separately validated by comparison with the few existing 
measurements |[30ll3Tl[32ll33ll34Tl .The main systematic uncertainty on the ITS standalone analysis comes 
from the tracking, due to the high occupancy and small number of tracking points. This was estimated 
from data using global tracks as a reference. The other systematic uncertainties are smaller, and include 
the effect of the magnetic field configuration (ExB effect), of the track selection and of the PID cuts. 
Similarly, the uncertainties related to the tracking efficiency in the TPC were investigated comparing 
different sets of tracks in data and Monte Carlo and by a variation of the track cuts. The uncertainty 
related to the combined TPC/TOF PID procedure was estimated varying the PID cut between 2 and 
5 a. The tracks reaching the TOF detector have to cross a substantial amount of additional material 
(X/X « 0.23), mostly due to the TRD US. Since the TRD was not fully installed in 2010, the analysis 
was repeated for regions in azimuth with and without installed TRD modules, allowing one to determine 
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Table 1: Main sources of systematic uncertainty. See text for details. 
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0.1 3 


0.2 3 


0.35 
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3% negl. 
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0.65 
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10% 


10% 




10% 


ITS PID 


2% 


4% 
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Pt range (GeV/c) 


0.3 0.65 


0.3 0.6 
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0.8 


plnHal tracking 
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4% 


4% 




4% 


TPC PID 


3% 


5% 




1.5% 


pi range (GeV/c) 


0.5 3 


0.5 3 


0.5 


4.5 


TOF matching 
efficiency 


3% 


6% 




3% 


TOF PID 


2% 7% 


3% 15% 


5% 


25% 



the uncertainty due to the aditional material. The systematics related to the PID extraction in the TOF 
analysis were estimated varying the parameters of the expected sources by ±10%. 

The pi distributions of positive and negative particles were found to be compatible within errors, we 
therefore show results for summed charge states in Fig. [T] The spectra are compared to RHIC results 
in Au-Au collisions at i/5nn = 200 GeV ll35l l36l and to hydrodynamic models. The spectral shapes 
show a significant change from RHIC to LHC energies, having a distinctly harder distribution. Within 
hydrodynamic models, this indicates a significantly stronger radial flow. In the range pi < 1.5 GeV/c 
VISH2+1 [37], a viscous hydrodynamic model, reproduces fairly well the pion and kaon distributions, 
but misses the protons, both in shape and absolute abundance. In this model, the particle yields are taken 
to be thermal at T c h = 165 MeV (see below). The difference is possibly due to the lack of an explicit 
description of the hadronic phase in the model. This interpretation is supported by the comparison with 
HKM j38l l39l . a similar model in which, after the hydrodynamic phase, particles are injected into a 
hadronic cascade model (UrQMD H01 |4T1 ). which further transports them until final decoupling. The 
hadronic phase builds additional radial flow, mostly due to elastic interactions, and affects particle ra- 
tios due to inelastic interactions. HKM yields a better description of the data. At the LHC, hadronic 
final state interactions, and in particular antibaryon-baryon annihilation, may therefore be an impor- 
tant ingredient for the description of particle yields |j42l |39l , contradicting the scenario of negligible 
abundance-changing processes in the hadronic phase. The third model shown in Fig.[T](Krak6w P31I441 ) 
introduces non-equilibrium corrections due to viscosity at the transition from the hydrodynamic descrip- 
tion to particles, which change the effective T c h, leading to a good agreement with the data. In the region 
pj < 3 GeV/c (Krakow) and pj < 1.5 GeV/c (HKM) the last two models reproduce the experimental 
data within ~20%, supporting a hydrodynamic interpretation of the transverse momentum spectra at the 
LHC. These models also describe correctly other features of the space-time evolution of the system, as 
measured by ALICE with charged pion correlations RBI . In order to quantify the kinetic freeze-out pa- 
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Fig. 1: (color online) Transverse momentum distributions of the sum of positive and negative particles (box: 
systematic errors; statistical errors smaller than the symbol for most data points), fitted individually with a blast 
wave function, compared to RHIC data and hydrodynamic models. 



rameters at ^Ann = 2-76 TeV, we performed a combined fit with a blast wave function [15]. It should 
be noted that the value of the 7jy„ parameter extracted from the fit is sensitive to the fit range used for 
the pions, because of the large contribution from resonance decays (mostly at low pj), which tend to 
reduce T^. For this reason, the pi ranges 0.5-1 GeV/c, 0.2-1.5 GeV/c, 0.3-3 GeV/c for %, K, and p 
were used. These hydro-motivated fits do not replace a full hydrodynamic calculation, but allow one 
to compare with a few parameters the measurements of different experiments. The data are well de- 
scribed by the combined blast wave fit with a collective radial flow velocity (j8t) = 0.65 ±0.02, and a 
kinetic freeze-out temperature of 7jy„ = 96 ± 10 MeV. As compared to fits to central Au-Au collisions 
at t/snn = 200 GeV/c, in similar pi ranges (35l @6l, (jSt) at the LHC is ~10% higher while Tu n is 
comparable within errors. 

The mid-rapidity < 0.5) pi -integrated particle yields were extracted by fitting the %, K, and p spectra 
individually with a blast wave function, in order to extrapolate to zero pj. The individual fits are shown 
in Fig.[T]as solid curves; the fraction of extrapolated yield is small: about 7%, 6%, and 4% for n, K, and 
p. Its uncertainty was estimated using different fit functions ll24l . The particle ratios are compared in 
Fig.|2]to results at ^/snn = 200 GeV and to the predictions from thermal models, using Hb = 1 MeV and a 
T c h of 164 MeV or 170 MeV [17]. The value for Hb is based on extrapolation from lower energy data. 
T c h was found to be constant above a center-of-mass energy of a few ten GeV, so the value obtained from 
fits to RHIC data was used. The systematic uncertainties on the particle ratios were computed taking into 
account the correlated sources of uncertainty (mainly due to the tracking efficiency for different particles 
and to PID and extrapolation for anti-particle over particle ratios). In the following we quote the total er- 
ror for the ratios, as the statistical error is negligible. The anti-particle/particle ratios are all unity within 
errors, consistent with a vanishing baryochemical potential In order to minimize the sensitivity to 
Hb, the ratios K/tt = (K + + K~)/ (ti + + 7i~) and p/n = (p +p)/ (7T + + 7i~) are also shown. The ratio 
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Fig. 2: (color online) Mid-rapidity particle ratios, compared to RHIC results and predictions from thermal models 
for central Pb-Pb collisions at the LHC (combined statistical and systematic errors). 



K/7T = 0.149 ±0.010, is similar to the lower energy values and agrees with the expectations from the 
thermal model [7]. However, the ratio p/n = 0.046 ±0.003, is significantly lower than expected, by a 
factor ~ 1.5-1. 9 (p/7T~ 0.07 — 0.09 for Q and lfl~71 respectively). The two models differ mainly in the 
hadron mass spectrum implementation, but were both successful in describing RHIC data. The compar- 
ison with RHIC data also hints at a slight decrease of the p/n ratio with energy (by a factor ~1.2), while 
essentially no change was predicted. The thermal models proved to be very successful over a wide range 
of energies (from ^/^ NN = 2 GeV to ^Ann = 200 GeV lHl|2l|6l[T0|): such a large difference for one of the 
most abundantly produced particle species was therefore unexpected. In retrospect, some disagreement 
between data and the thermal model is also apparent in the RHIC data, with the proton measurements 
being about 20% lower than predictions flU @6l @71. However, this difference was not considered to 
be significant, because of the differences between model implementations, model uncertainties [48 ] and 
experimental uncertainties in the subtraction of secondary particles in the RHIC experiments. This is- 
sue will likely be clarified by a thermal analysis including strange and multi-strange baryons at the LHC. 
Current speculations are that final state interactions in the hadronic phase, in particular via the large cross 
section channel for antibaryon-baryon annihilation 11421 . could explain the significant deviation from the 
usual thermal ratios. A similar conclusion is implied by the HKM model, where p/n = 0.052, consistent 
with our measurement [39]. An alternative scenario conjectures the existence of flavor and mass depen- 
dent pre-hadronic bound states in the QGP phase, as suggested by recent lattice QCD calculation and 
QCD-inspired models Il49ll50l. 

In summary, we presented the first measurements of 71, K, and p production in central Pb-Pb collisions at 
Y^nn = 2.76 TeV at the LHC. The pj distributions are harder than previously measured at RHIC. They 
are well described by hydrodynamic models including a refined description of the late fireball stages. 
Fitting the spectra with a hydro-inspired blast wave model results in the highest radial flow parameter 
ever measured, (j&r) =0.65 ± 0.02. The integrated particle ratios were compared with expectations from 
thermal models. While the K/n ratio was found to agree with these expectations, p/n is a factor > 1.5 
lower. 
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